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I.  INTRODUCTION 

Since  one  can  quite  easily  and  accurately  compute  the  rotational  and  vibrational-rotational 
absorption  spectra  of  diatomic  molecules  directly  from  quantum  mechanics,  it  would  at  first 
glance  seem  useless,  even  quixotic,  to  try  to  approximate  these  spectra  by  means  of  a  much 
more  laborious  classical  technique.  The  justification  for  such  an  exercise  is  that  this  classical 
approach  may  be  extended  to  systems  in  which  a  quantum  approach  is  as  yet  impractical.  For 
example,  one  may  compute  the  infrared  spectra  for  systems  involving  many  atoms,  such  as 
large  molecules,  clusters  or  liquids,  in  which  a  normal  mode  analysis  may  not  be  feasible.1  In 
addition,  one  can  extend  the  approach  to  explicitly  time  dependent  problems  such  as  the  com¬ 
putation  of  the  transient  vibrational  spectrg  one  might  hope  to  measure  during  the  course  of  a 
chemical  reaction.1'3  It  is  thus  useful  to  test  this  approach  with  a  simple  molecular  system  for 
which  the  results  can  be  rigorously  tested  both  against  accurate  quantum  calculations  and 
against  experimental  measurements. 

In  Fig.  1  we  see  the  rotational  and  fundamental  vibrational-rotational  absorption  bands  of 
CO  which  we  will  use  for  our  test  of  the  ability  of  classical  mechanics  to  reproduce  quantum 
reality  in  the  form  of  rotational  and  vibrational  spectra.  The  shapes  of  vibrational  spectra  are  a 
sensitive  reflection  of  both  inter-  and  intra-moiecular  motion4  and  can  thus  in  principle  supply  a 
method  for  discovering  the  microscopic  trajectories  responsible  for  chemical  processes  in  solu¬ 
tion. 

II.  THEORETICAL  TECHNIQUES 

The  theoretical  approach  consists  of  five  steps.1  First,  classical  mechanics  is  used  to  com¬ 
pute  the  atomic  trajectories  for  the  system  of  interest  from  a  set  of  initial  atomic  coordinates 
and  momenta  and  a  given  potential  surface.  Second,  a  time-varying  dipole  moment  vector  (for 
infrared  spectra)  or  a  time-varying  polarizability  matrix  (for  Raman  spectra)  for  the  entire  sys¬ 
tem  is  calculated  from  the  trajectories.  Third,  linear  response  theory4*6  is  used  to  derive  an 
infrared  (or  Raman)  spectrum  "specific"  to  the  chosen  initial  coordinates  and  momenta  from 
the  power  spectrum  of  the  dipole  moment  (or  polarizability)  time  histories.  Fourth,  the 
"specific"  spectra  are  averaged  over  the  ensemble  of  initial  coordinates  and  momenta  which  is 
appropriate  to  the  experimental  conditions  of  interest  (e.g.  a  particular  temperature  and  pres¬ 
sure).  Fifth,  quantum  corrections  are  applied  where  necessary  to  converge  toward  quantum 
reality. 

By  these  five  steps,  surprisingly  accurate  infrared  and  Raman  spectra  may  be  derived  from 
three  functions  of  atomic  coordinates:  potential  energy,  dipole  moment  vector  and  polarizabil¬ 
ity  matrix.  We  will  illustrate  here  the  electric  dipole  (microwave  and  infrared  absorption  and 
emission)  case,  and  leave  the  details  of  Raman  scattering  for  another  paper. 
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A.  Molecular  Dynamics  and  Classical  Trajectories 

Our  first  step  is  to  choose  a  set  of  initial  coordinates  and  momenta  Tor  the  collection  of  N 
atoms  from  an  ensemble  of  possibilities  representative  of  the  system  of  interest,  for  example  a 
system  at  a  particular  pressure  and  temperature.  We  make  this  choice  purely  classically 
although  Marcus  and  co-workers  have  shown7  how  one  might  also  make  such  choices  semiclas- 
sically.  The  second  step  is  to  compute  the  classical  trajectories  rt (/), . . . ,  rA(/)  describing  the 
atomic  motions  by  integrating  Newton’s  Second  Law 
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in  which  V  -  V(rl, . . .  »r,v)  is  the  potential  energy  of  the  atoms  at  positions 
r|f . . .  ,r,v,  F,  -  F,(r|, . . .  ,r,v)  is  the  force  on  the  /  th  atom,  and  m,  is  the  mass  of  the  i  th 
atom.  A  modified  Verlet  integration  algorithm  is  used.8'9  For  the  example  shown  below  of  a 
CO  solution  in  Ar,  minimum  image  periodic  boundary  conditions  with  truncated  octahedral 
boundaries10  are  used  to  reduce  edge  effects. 


B.  Time  History  of  Dipole  Moment 

Given  the  function  n(rlt . . . ,  r,v) ,  which  describes  the  dipole  moment  of  the  system  as  a 

function  of  the  atomic  positions,  we  can  use  the  trajectories  of  the  atoms,  r,(/) . r*(r),  to 

compute  /«.(/),  the  dipole  moment  of  the  system  as  a  function  of  time. 


C.  Linear  Response  and  Specific  Spectra 

For  a  system  at  equilibrium,  we  can  use  linear  response  theory4-6  to  compute  the  infrared 
spectrum  o(o»)  from  the  dipole  moment  time  history,  /t(f).  We  can  thus  relate  the  spectrum 
of  the  fluctuations  which  /i(f)  naturally  undergoes  at  equilibrium  to  the  spectrum  with  which 
jt (f)  responds  when  driven  by  an  external  oscillating  electric  field,  i.e.  we  relate  the  equilibrium 
dipole  moment  fluctuation  spectrum  to  the  absorption  spectrum  when  irradiated  with  light. 

The  appropriate  linear  response  equations  are4-6 

«(*,)  -  /(&>)  (2) 
iTicn 

09 

/(«)  =  (2rr)_l  J dt  exp(— i<ot)  <§i(0)-/t(t)>  (3) 


in  which  a( w)  is  the  absorption  cross  section  as  a  function  of  angular  frequency  w,  £  = 
(kBT)~l  in  which  kB  is  Boltzmann’s  constant  and  T  the  temperature,  1i  =  h/2iz  in  which  h  is 
Planck’s  constant,  c  is  the  speed  of  light,  n  is  the  index  of  refraction  of  the  medium,  /(<u)  is 
defined  as  the  absorption  lineshape  and  is  evaluated  for  an  isotropic  medium,  and 
<ft(0)-p(/)>  is  the  ensemble  average  of  the  dipole  moment  time  correlation  function. 


As  Parseval’s  theorem11  and  the  Wiener- Khintchine  theorem  6-12  show,  /(«>)  may  be 
computed  in  the  mathematically  equivalent  form  of  the  power  spectrum7 


/  (o»)  -  (2jt)  1 


limy-  £  exp(-/WV,(f)lJ 
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allowing  the  use  of  fast  Fourier  techniques.  Considerable  care  must  be  used  to  properly  apply 
spectral  estimation  theory,  windowing  and  windowing  corrections,11- 12  in  both  the  correlation 
and  power  spectra  methods  to  avoid  distorting  the  spectra  as  a  result  of  the  use  of  finite  time 
histories.  Particular  care  must  be  taken  if  band  wings  are  to  be  correct  or  if  small  bands  are  to 
be  seen  in  the  presence  of  large  ones.  We  use  a  four  term  -74  db  Blackman-llarris  window,1-1 
which  allows  the  correction  to  be  applied  as  a  simple  frequency-space  convolution  of  the 
Fourier  transform  of  the  dipole  moment  time  histories  with  the  Fourier  transform  of  ihc.win- 
dow  function.  This  is  done  before  computing  the  power  spectrum.  The  windowed  power 
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spectrum  is  multiplied  by  the  inverse  of  the  sum  of  the  squares  of  the  windowing  function 
which,  by  Parseval's  theorem,11  can  be  evaluated  either  in  time  (windowing  function)  or  fre¬ 
quency  (Fourier  transform  of  windowing  function).  This  simple  scaling  factor  corrects  the 
spectral  band  areas  for  the  scaling  effect  of  the  windowing  and  is  correct  because  of  the  loss  of 
phase  coherence  due  to  the  ensemble  averaging.  The  windowing  and  the  area  scaling  correction 
can  be  combined  into  a  convolution  of  the  Fourier  transform  of  the  dipole  moment  time  his¬ 
tory  with  the  discrete  function  fk  =»  (-0.00170,  0.08725,  —0.46174,  0.74724,  —0.46174, 

0.08725,  -0.00170)  for  k  -  -3 . 0, . . .  ,+3  before  taking  the  square  of  the  absolute  value 

in  computing  the  power  spectrum. 

D.  Ensemble  Averaging 

The  specific  spectrum  computed  from  a  single  time  history  is  the  spectrum  of  a  particular 
very  small  sample  for  a  particular  very  short  time  period.  To  compare  to  an  experimental  spec¬ 
trum,  an  ensemble  average  must  be  made  over  a  distribution  of  specific  spectra  corresponding 
to  the  experimental  conditions.  For  example,  to  compare  with  a  spectrum  taken  at  temperature 

Tone  can  choose  a  reasonable  initial  guess  for  the  positions  rj C0> . r,v(0)  of  the  atoms, 

and  choose  the  velocities  ^(0),  — ,  f,v(0)  randomly  from  a  Maxwell-Boltzmann  distribution 

for  temperature  T.  After  computing  the  trajectories  rt (/) . r.v(r)  forward  in  time  for  a 

.  suitable  period  to  equilibrate  the  system  (in  particular  to  equilibrate  potential  and  kinetic  ener¬ 
gies)  one  can  preserve  positions  but  change  to  a  new  set  of  velocities,  again  chosen  randomly 
from  the  Maxwell-Boltzmann  distribution,  and  then  continue  integrating  forward  in  time. 
(Such  velocity  randomization  is  important  to  achieve  a  proper  ensemble  average  for  low  density 
gas  phase  molecules  for  which  energy  and  angular  momentum  are  otherwise  conserved  over 
long  periods.)  Several  repetitions  of  this  process  will  result  in  a  system  equilibrated  at  the 
desired  temperature  T  and  randomly  chosen  as  to  velocity.  The  positions  of  the  atoms  may 
still,  however,  be  correlated  with  the  initial  guessed  positions.  To  sample  more  of  configuration 
space  we  therefore  use  the  final  positions  of  each  run  as  the  initial  guess  for  the  next  run  and 
again  repetitively  randomize  velocities.  An  alternative  method  would  be  to  use  the  semidassi- 
cal  techniques  of  Marcus  and  co-workers7  and  quantize  the  action  variables  and  then  perform 
the  ensemble  average  over  a  Boltzmann  population  of  these  quantized  systems. 

The  spectra  which  are  derived  from  classical  molecular  dynamics,  classical  linear  response 
theory,  and  classical  statistical  mechanical  ensemble  averaging  we  will  call  "Newtonian."  They 
are  straightforward  to  calculate  in  principle,  but  can  be  unusually  demanding  arithmetically. 
Thus  for  their  computation  we  use  what  might  be  called  an  "instrument  for  theory,"  a  network 
of  processors  (including  an  array  processor)  and  generalized  program  package,  which  are 
described  elsewhere.1-2 

E.  Quantum  Corrections 

In  order  to  discover  what  corrections  are  necessary  to  bring  our  classical  calculations  into 
agreement  with  quantum  reality,  we  will  compare  three  different  ways  of  computing  the  absorp¬ 
tion  spectrum:  /)  quantum  mechanics  without  linear  response  theory,  ii)  correspondence  princi¬ 
ple  classical  mechanics  without  linear  response  theory  as  given  by  the  limit  of  quantum 
mechanics  as  Planck's  constant  approaches  zero,  and  Hi)  ordinary  Newtonian  classical  mechan¬ 
ics  with  classical  linear  response  theory  and  ensemble  averaging  from  classical  statistical 
mechanics. 
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1.  Quantum  solution 

Quantum  mechanically,  we  can  evaluate  the  absorption  spectrum,  aim),  from  l  (at)  as 
shown  in  Eq.  (2),  without  recourse  to  linear  response  theory  as,4-6 

/(at)  -  X  Z  Pi  l</l#»l'>l2  S(wy,  -  (i>)  (5) 

i  l 

in  which  we  have  averaged  over  an  isotropic  distribution  of  molecular  orientations.  In  the 
above,  i  is  an  initial  quantum  state,  /  a  final  state,  p,  -  exp(-/3£)/]Texp(-/3£,)  is  the  proba- 

I 

bility  for  initial  state  /,  </lp|/>  is  the  electric  dipole  transition  matrix  element  between  states 
j  and  /,  a>),  -  (Ef—  E,)/1I,  E,  and  £,  are  the  energies  of  the  final  and  initial  states,  respec¬ 
tively,  and  $  is  the  Dirac  delta  function.  In  the  Appendix,  we  give  the  evaluation  of  this 
expression  explicitly  for  the  case  of  the  isolated  diatomic  molecule  and  thus  arrive  at  a  quantum 
mechanical  computation  of  the  rotational  and  vibrational-rotational  spectrum  a(at),  by  inserting 
the  resulting  /(at)  into  Eq.  (2). 


2.  Correspondence  principle  classical  solution 

We  can  solve  problems  in  classical  mechanics  in  either  of  two  very  different  ways.  We 
can  take  the  usual  route  of  Newton’s  laws,  or,  as  in  this  section,  we  can  solve  the  problem 
quantum  mechanically,  but  take  the  correspondence  principle  limit  as  7T — 0.  Both  forms  of 
classical  mechanics  are  useful  in  bringing  different  viewpoints  to  bear  on  understanding  the 
problem,  and  in  providing  checks  of  theoretical  and  computational  methods.  The  correspon¬ 
dence  principle  approach,  i.e.  the  limit  of  infinitesimal  photons  causing  infinitesimal  perturba¬ 
tions  of  the  molecular  system,  gives  for  the  classical  7f— 0  limit  of  Eq.  (2) 

a(at)  .  i2d2f!£  /(<a)  (6) 

3c/i 

which  we  use  with  our  classical  molecular  dynamics  calculations. 

However,  in  the  correspondence  principle  calculations  illustrated  in  the  figures,  we  have 
taken  the  /i— 0  limit  by  iterating  each  calculation  of  a  (at)  using  successively  lower  values  of  7F 
in  Eqs.  (A7-A9)  of  the  Appendix  and  in  Eq.  (2)  until  the  computed  a  (to)  converges,  which 
occurs  by  the  time  U  has  been  replaced  by  —0.0177. 


3.  Evaluation  of  quantum  corrections 

Simple  quantum  corrections  are  used  to  bring  the  classically  computed  spectra  into  agree¬ 
ment  with  theoretical  and  experimental  quantum  reality.  The  correction  to  the  rotational  band 
contour  will  be  seen  to  be  minor  and  usually  negligible.  The  vibrational-rotational  contour 
correction,  while  not  negligible,  is  relatively  small  and  easy  to  apply,  consisting  of  a  frequency 
shift  and  a  shape  correction. 

We  write  Eq.  (2)  as 

«<?(«)  -  ,Q{u)  (7) 

377  ch 

in  which  the  superscript  Q  indicates  quantum.  If  we  substitute  -at  for  at,  interchange  the  sub¬ 
scripts  i  and  /  in  Eq.  (5),  and  recognize  that  |</||tl'>l2  "  I</ImI/>Ij  and  S(w)  -  S(-«), 
we  then  find  that 

/Q(-at)  —  £  iLPy  I  </\n\i>  lJ  hCatji—at).  (8) 

(  I 

For  a  system  at  equilibrium,  p,  -  exp(07Ta >„)pi,  and  thus 

lQ(at) 
l°(-at) 


exp(/J7rw„)  -  exp(/37Ia>)  , 


(9) 
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which  is  just  an  expression  of  detailed  balance.  Because  of  the  8-functions  in  Eqs.  (S)  and  (8), 
<i*  can  replace  <u/,.  The  correspondence  principle  classical  limit  is  thus 

Ic(m)  _  Mu,) 


lc(-u>) 


lim 
w — o 


-  1 


(10) 


in  which  the  superscript  C  indicates  classical.  Eqs.  (9)  and  (10)  can  be  understood  as  follows. 
The  energy  carried  by  a  photon  of  angular  frequency  u>  connects  states  in  a  quantum  system 
which  are  Hut  apart  and  of  appreciably  different  Boltzmann  probabilities  in  an  equilibrium  sys¬ 
tem,  and  thus  the  /tf(w)  for  absorption  and  l°(-ut)  for  emission  are  related  by  the  detailed  bal¬ 
ance  factor  exp(/37Ft «).  In  the  classical  correspondence  limit  as  7/ — 0,  the  photon  connects 
states  of  only  infinitesimally  different  energy  thus  lc{ut)  for  absorption  and  /c(-a*)for  emission 
become  equal.  Comparison  of  Eqs.  (9)  and  (10)  suggests  that  one  should  consider  the  quan¬ 
tum  correction5- 14<  15 


/c(l!)  ”exp(^w/2) 


(11) 


which  symmetrically  and  simply  fulfills  the  requirements  of  both  equations.  The  factor  of  ut/2 
inside  the  exponential  in  Eq.  (11)  arises  because  lQ(-ut)  and  lQ(at),  which  are  related  by 
detailed  balance,  lie  2w  apart  in  angular  frequency. 

a.  Rotational  correction.  The  classical  absorption  cross  section  ac(w)  is,  from  Eq.  (7), 


ac(ut)  “  lim  a^(ut) 
»— o 
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[1— exp(— fffffai)] 
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-  Mut).  (13) 

3cn 

in  which  we  define  the  correspondence  principle,  7j-*0,  classical  limit  of  lQ(ut)  as  lc(u>). 

Using  Eqs.  (7),  (11)  and  (13)  we  have, 

gQ(ut)  sinh(fflTm/2) 

ac(o.)  “  pHut/2 

which  can  be  used  as  a  quantum  correction  to  ac(ut),  the  classical  pure  rotational  band  contour. 
Since  for  small  jc,  sinhU)  =  x,  the  quantum  correction  is  approximately  unity  for  low  frequen¬ 
cies  and  high  temperatures,  and  is  only  1.0096  at  100  cm-1  and  300  °K.  Thus,  in  general  its 
effect  on  the  rotational  spectrum  is  small  and  the  quantum  and  classical  pure  rotational  band 
contours  are  essentially  the  same. 

b.  Vibrational-rotational  correction.  To  derive  quantum  corrections  for  the  vibrational- 
rotational  spectra,  we  first  make  the  approximation  of  separating  vibration  and  rotation  (a 
separation  not  made  in  our  molecular  dynamics  which  can  mix  and  couple  translational,  rota¬ 
tional  and  vibrational  motions),  giving  for  the  matrix  element  in  Eq.  (5),  as  shown  in  Eq.  (A4) 
of  the  Appendix, 

</ImI'>  “  <v'J'M'\n\vJM>  (15) 

=  <v’|m(K)|v>  <J’M'\R\JM>  (16) 

in  which  n(R)  is  the  magnitude  of  the  dipole  moment  and  R  is  a  unit  vector  along  the  inlemu- 

dear  axis.  We  also  drop,  for  example  in  Eq.  (All),  the  coupling  terms  between  vibration  and 
rotation  in  the  energy  which  gives 


•Hu>fi  -  f(v'.y')  -  E(v,J) 

=  !£(v')  -  £(v)]  +  [£U')  -  £(/))  —  Tfw, ,  +  Hutjj  . 


(17) 


where  «vr  and  utjj  are  angular  frequency  differences  between  states  »•'  and  v  and  J'  and  J 
respectively.  Substituting  Eqs.  (16)  and  (17)  into  Eq.  (5)  we  get 


(18) 


/°(«)  -  I  III  pM  PiU)  I</ImO?)|w>  <J'M'\R\JM>\2 

*  *•  J.MJM 


X  8  (<■»„'„  +  01  jj  —  tu). 


in  which,  for  a  system  at  equilibrium, 
.  _  exp[-gf(v)l 


p,(v)  - 


Iexpl-)8£(v)) 


„  , ,)  .  txpl-HE(J)] 
plU  lexpl -pE(J)\  • 

i 

Separation  of  the  sums  over  the  vibrational  and  rotational  states  gives 
lQ(a>)  -  [llp,(v)  |<vV(*)|v>|2|  [l  I  P,U)  \<J'M'\R\JM>\‘ 

I  *  V  I  \j.MJM 

x  S(oivv  +  oi jj  —  oi)  . 


We  can  write  this  in  terms  of  a  convolution16 
/°U)  -  III  p,(v)  |<v'|M(/?)|v>|J  8(«vV  -  a,) 


*  I  I  PiU)  \  <JM'\R\JM>\2  8 (o,jj  -  <o) 

J.M  J.M 


-  I?(oi)  *  /£(«), 


(19) 

(20) 

(21) 

(22) 

(23) 


in  which  the  quantum  absorption  lineshape  /°(w)  is  seen  to  be  a  convolution  of  the  quantum 
vibrational  absorption  lineshape  l?(oi)  and  the  quantum  rotational  absorption  lineshape  /«(w), 
if  vibration  and  rotation  are  separable  by  Eqs.  (16)  and  (17).  Therefore,  in  this  approximation 
of  separation  of  vibration  and  rotation,  we  can  quantum  correct  the  vibrational  and  rotational 
lineshapes  f?(oi)  and  /£(«)  separately  and  then  convolve  the  results. 

Similarly,  when  vibration  and  rotation  are  separable  by  Eqs.  (16)  and  (17),  by  applying 
the  frequency  convolution  theorem11  to  Eq.  (23),  we  can  separate  the  dipole  moment  time 
correlation  function4-6  into  a  product  of  vibrational  and  rotational  correlation  functions17 


<#i(0)jt(r)>  -  J* dm  exp (iuit)  /<>(«)  «*  J* du  exp(fW)[  l?(oi)  *  /£(<u)| 

—  SO 

-  Jjrfw  expO'ca/)  /?(«)]  |  J  du  exp(/wr)/f(w) 


-  <M(0 )pU)>  <R(0)R(r)> 


Substituting  Eq.  (23)  into  Eq.  (7),  the  quantum  vibrational-rotational  absorption  cross 
section,  i.e.  spectrum,  is 


and  the  classical  spectrum  is 
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ac(«)  -  lim  a^(  Itf) 

7f—-0 


4it2ci)  ||  1— exp(— 07}<o) 

3  cn  n — o  Tf 


lim  /M  *[lim  /«(«<) 
,n-o  [n-o 


(28) 


4n2(u2(3 

3cn 


(/£(«)  *  /#(oj)]  . 


(29) 


A  comparison  of  Eqs.  (27)  and  (29)  shows  that  given  the  classical  vibrational-rotational 
spectrum,  ac(<o),  and  assuming  separability  of  vibration  and  rotation,  we  can  quantum  correct 
it  to  approximate  the  quantum  band  contour  by  the  following  three  steps. 

i)  Initial  factor.  The  ratio  of  the  quantum  and  classical  versions  of  the  first  factor  in 
brackets  in  Eqs.  (27)  and  (29)  is 


4 ir2<u  _  [1 — exp( — g7T cu)] 

3cn _ 5 _ _  1— exp(— ffftt») 

4ir2u20  ptifu 

3cn 


(30) 


and  we  multiply  the  classical  spectrum  by  this  factor. 

ii)  Vibrational  linesliape.  From  Eqs.  (22)  and  (23),  we  have 
/?(«)- 15>,(v)  |<  v'|/tt(/?)|v>|2  5(a»„„  —  <u)  .  (31) 

V  V* 

If  we  evaluate  /P(w),  the  vibrational  absorption  lineshape  function,  using  the  linear  dipole 
function  approximation,  p.(R)  =  /a0  +  Pi^R.  for  an  equilibrium  system  of  harmonic  oscillators 
(see  Appendix)  we  find18 


/?(*>) 


_ _ 

2mrw[\—  exp(-/J7T<u)] 


M?  8(wff-iu) 


(32) 


in  which  mr  is  the  reduced  mass  of  the  oscillator,  <i>„-(A/m,)'*  and  k  is  the  force  constant. 
Taking  the  correspondence  principle  limit  of  Eq.  (32)  yields 


If  (.to)  ■*  lim  /?(„)  —  (2wr/3<o:)  1  fj.}  &(<i>0— co)  , 

Jfl — 0 

giving  a  quantum  correction  factor  from  Eqs.  (32)  and  (33)  of 
l£(<o)  1-exp 


(33) 

(34) 


This  harmonic  oscillator  quantum  correction  factor  exactly  cancels  out  the  quantum  correction 
in  Eq.  (30)  above.  Therefore,  for  a  harmonic  oscillator  with  linear  dipole  moment  function, 
there  is  no  quantum  correction  and  the  quantum  and  classical  spectra  agree  for  pure  vibration. 
The  quantum  corrections  to  vibration  are  thus  due  to  anharmonicity  in  the  potential  and  non¬ 
linearity  in  the  dipole  moment  function.  We  will  discuss  only  the  potential  anharmonicity  here. 

For  the  fundamental  absorption  band,  the  sum  of  all  v+1— v  transitions,  the  dominant 
effect  of  potential  anharmonicity  for  pure  vibration  is  a  shift  in  the  band  center  to  a  lower  fre¬ 
quency.  That  such  a  quantum  correction  is  needed  can  be  seen  qualitatively  as  follows.  For  a 
usual  high  frequency  vibration  at  room  temperature,  the  classical  oscillator  is  at  —  kBT,  near  the 
bottom  of  the  well  in  a  largely  harmonic  region,  while  the  quantum  oscillator  is  transitioning 
largely  from  the  zero  point  v=0  to  v-1  levels  (both  at  several  kHT)  thus  feeling  much  more  of 
the  anharmonicity. 

A  simple  correction  suffices  for  the  Av-1  fundamental  transition  treated  here.  We  need 
only  shift  the  whole  band  as  calculated  classically  over  to  its  proper  quantum  location,  without 
needing  to  alter  its  shape.  This  can  be  done  by  computing  the  shift  of  the  rotaiionlcss  oscillator, 
involving  the  fictitious  {\\J)  transitions  (v+1,0)— (r.O).  We  can  derive  this  pure  vibrational 
spectrum  from  Eq.  (A8)  in  the  Appendix  as 
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/?(w)  -  I  p(v, 0)  p?  l(v+l)/2a]  s|-^±-)’^-.£(^0)  -  J  .  (35) 

,mO  n 

We  evaluate  the  center  in  angular  frequency  of  aim)  for  the  vibrational  band  as  computed  from 
Eqs.  (2)  and  (35),  first  for  the  quantum  case  and  then  for  the  classical  case  by  reducing 
Planck’s  constant  until  convergence  is  reached.  This  offset  between  the  classical  and  quantum 
rotationless  band  centers  is  then  used  to  shift  the  classical  vibrational-rotational  band  contour 
into  its  quantum  corrected  position. 

While  this  procedure  can  be  quickly  carried  out,  there  is  also  a  rougher  approximation 
which  can  be  used.  Classical  oscillation  is  often  so  close  to  the  bottom  of  the  potential  well  that 
a  classical  vibration  sees  essentially  only  the  harmonic  part  of  the  potential  and  thus  only  the 
term  in  Eq.  (All)  in  the  Appendix  is  important.  Quantum  mechanically,  if  /?»>,.  is  high  enough 
with  respect  to  kBT,  only  v-0  will  be  appreciably  populated,  and  only  1 — 0  vibrational  transi¬ 
tions  need  be  considered.  Then  the  offset  in  angular  frequency  m  is  approximately  -4iri>exe,  as 
may  be  computed  from  Eq.  (Al  1). 

iii)  Rotational  lineshape.  Comparing  Eqs.  (27)  and  (29),  and  noting  Eq.  (11),  we  can 
quantum  correct  the  rotational  absorption  lineshape  function  /«( w)  for  the  vibrational- 
rotational  band  by  multiplying  it  by  exp(/37T(Aw)/2),  in  which  Aw  is  the  frequency  measured 
from  the  rotationless  center  of  the  vibrational  band.19  -0  This  factor  can  be  understood  by  con¬ 
sidering  rotation  as  a  mechanically  and  statistically  separate  process  from  vibration  in  which 
+Aw  is  a  rotational  absorption  from  /*— /  and  —Aw  a  rotational  emission  from  f—i,  and  the  fre¬ 
quency  difference  2Aw  therefore  separates  processes  for  which  detailed  balance  gives  a  probabil¬ 
ity  ratio  of  exp[07i(Aw)]  at  equilibrium. 

Thus,  the  classical  fundamental  vibration-rotation  spectrum  can  be  quantum  corrected  by 
two  easy  steps.  First,  the  spectrum  is  frequency  shifted  to  take  into  account  the  effect  of  vibra¬ 
tional  anharmonicity,  and  second,  it  is  multiplied  by  exp[/J7/(Aw)/2].  This  assumes  that  vibra¬ 
tion  and  rotation  are  approximately  independent,  so  that  detailed  balance  may  be  applied 
separately  to  each.  This  assumption  is  not  necessarily  appropriate  for  condensed  systems. 

In  Fig.  2  we  show  some  of  the  effects  of  the  various  quantum  corrections  discussed  above 
for  the  fundamental  vibration-rotation  absorption  band  of  gas  phase  CO.  First  of  all,  the  upper 
panel  shows  the  quantum  band  contour  from  Eqs.  (2),  (5)  and  (A8-A18)  as  a  solid  line.  The 
dotted  line  shows  the  classical  band  contour  from  the  7/— 0  limit  of  these  equations.  The 
dashed  line  shows  the  effect  of  applying  the  vibrational  anharmonic  offset  quantum  correction 
of  -22  cm-1  at  300 TC  from  Eq.  (35)  (the  rougher  approximation  of  —4tti>cxc  gives21 
—26  cm-1).  The  dashed  and  dotted  line  shows  the  effect  of  applying  the  rotational  detailed  bal¬ 
ance  correction  exp[(J7F(Aw)/2l  alone  to  the  classical  vibration-rotation  band  contour.  The 
effect  plying  both  the  vibrational  offset  and  the  rotational  detailed  balance  corrections  is 
shown  •  .pen  circles,  and  is  seen  to  match  very  closely  the  true  quantum  contour. 

The  lower  panel  of  Fig.  2  shows  the  large  effect  of  vibrational  anharmonicity  on  the  band 
contour.  The  computations  are  similar  to  those  shown  in  the  upper  panel,  except  that  CO  is 
treated  as  a  harmonic  oscillator.  The  same  equilibrium  internuclear  distance  and  second  deriva¬ 
tive  of  the  potential  are  used  as  for  previous  case,  but  the  third,  fourth,  and  higher  derivatives 
are  set  to  zero.  There  is  now  no  anharmonic  offset.  The  solid  line  shows  the  quantum  band 
contour,  the  dotted  line  the  classical  contour,  and  the  open  circles  the  effect  of  applying  the 
rotational  detailed  balance  correction.  It  should  be  noted  for  the  real  anharmonic  CO  band  con¬ 
tour  that  the  difference  in  peak  intensity  of  the  P  and  R  branches  is  much  more  pronounced 
than  for  the  harmonic  case.  The  inability  to  match  the  P  and  R  branch  asymmetry  of  real  spec¬ 
tra  with  harmonic  oscillator  -  rigid  rotor  models  has  been  a  source  of  confusion.  It  appears  from 
Fig.  2  that  the  problem  has  not  been  a  quantum  effect  at  all,  but  simply  due  to  ignoring  the 
true  anharmonicity  of  the  potential. 


III.  PURE  ROTATIONAL  SPECTRA 

Fig.  3  shows  the  pure  rotational  band  for  gas  phase  CO.  The  upper  panel  shows  the  quan¬ 
tum  mechanical  spectrum,  computed  as  shown  in  the  Appendix  in  Eqs.  (A7)  and  (AI0-AI8) 
by  explicitly  calculating  the  individual  transitions  with  the  usual  time-dependent  perturbation 
theory  approach  to  compute  /(w),  and  then  using  Eq.  (2)  to  compute  «(<o).  The  solid  line 
shows  the  quantum  band  contour,  given  by  broadening  the  individual  transition  lineshapes  until 
they  merge,  using  Eq.  (A18).  The  CO  potential  used  in  all  these  calculations  is  Huffaker's 
analytical  expansion,-2- 25  giving,  for  use  in  the  quantum  and  correspondence  principle  classical 
calculations,  the  potential  derivative  values  at  the  equilibrium  internuclear  distance  R,.  of 
r(/?,)«1.92024xl0-3  J  nr2.  *""(/?,.)— 1.36445xl0“14  J  m“J  and  r"(Kr)-8.08245*  10“24 
J  m~4.  We  used21  1. 1 28323 x  10“ 10  m  and  reduced  mass  /nr«  1.1 3843x10“ 26  kg. 

The  middle  panel  of  Fig.  3  demonstrates  that  the  classical,  correspondence  principle,  h— *0 
limit  pure  rotational  band  contour  falls  almost  exactly  on  the  quantum  contour.  The  total  area 
of  the  rotational  band  is  1.88xi0“22m,  close  to  the  area  of  l.94xl0~::m  indicated  by  the  free 
rigid  rotor  sum  rule 

2  ir 

3xl07  '  /  ’ 

given  by  Buontempo  el  a/.24  Eq.  (36)  is  given  in  SI  units,  a-  is  the  absorption  cross  section,  X-1 
the  inverse  wavelength,  /a0  the  permanent  dipole  moment  and  /  the  moment  of  inertia. 

The  lower  panel  of  Fig.  3  demonstrates  that  "Newtonian"  classical  mechanics,  calculated 
from  classical  molecular  dynamics  as  in  Eq.  (1),  classical  linear  response  theory  as  in  Eqs.  (4) 
and  (6)  and  classical  statistical  mechanical  averaging  also  agrees  with  quantum  reality.  In  addi¬ 
tion,  it  is  seen  that  the  calculations  essentially  agree  with  the  experimental  gas  phase  measure¬ 
ments  of  Buontempo  et  a/24  on  1.6%  CO  solution  broadened  by  Ar  at  17.4  atm. 

Thus,  for  the  gas-phase  CO  rotational  band  contour,  /)  quantum,  //)  correspondence  prin¬ 
ciple  classical,  ///)  Newtonian  classical  and  h)  experiment  all  agree. 
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IV.  VIBRATIONAL-ROTATIONAL  SPECTRA 

Figure  4  shows  the  fundamental  vibrational-rotational  band  for  gas  phase  CO.  The  top 
panel  gives  the  results  from  quantum  mechanics,  the  total  effect  of  all  the  vibrational-rotational 
transitions  calculated  from  time-dependent  perturbation  theory  as  discussed  in  the  Appendix. 
The  quantum  vibrational-rotational  lines  in  the  P  and  R  branches  are  computed  explicitly  from 
the  quantum  transitions  and  then  checked  against  accurate  line  position  measurements.25  The 
quantum  band  contour  is  produced  by  widening  out  the  lines  as  gaussians  until  they  merge,  as 
calculated  from  Eq.  (2)  and  Eqs.  (A8-A18)  in  the  Appendix. 

The  middle  panel  of  Fig.  4  shows  that  the  correspondence  principle  classical  result,  as  cal¬ 
culated  as  the  7f— 0  limit  of  the  above  quantum  transition  result,  agrees  almost  exactly  with  the 
quantum  result  when  the  rotational  and  vibrational  quantum  corrections  described  above  are 
applied. 

The  lower  panel  of  Fig.  4  shows  the  Newtonian  classical  result,  computed  from  classical 
molecular  dynamics  by  Eq.  (1),  classical  linear  response  theory,  and  classical  ensemble  averag¬ 
ing,  and  then  quantum  corrected  in  the  same  way  as  in  the  middle  panel.  Newtonian  classical 
mechanics  is  seen  also  to  agree  very  closely  with  the  quantum  result.  In  addition,  as  shown, 
the  experimental  results  of  Armstrong  and  Welsh26  for  CO  broadened  by  104  amagats  of  He 
also  agree  in  shape  (no  experimental  absolute  intensity  was  given). 

Thus,  four  approaches  all  give  essentially  the  same  answer:  /')  quantum  mechanics  (time- 
dependent  perturbation  theory),  ii)  quantum-corrected  correspondence  principle  classical 
mechanics  as  the  1i—0  limit  of  quantum  mechanics.  Hi)  quantum-corrected  Newtonian  classical 
mechanics  plus  classical  linear  response  theory  and  classical  statistical  mechanics  and  /»•)  experi¬ 
mental  measurement. 


In  the  above  calculations,  we  used  a  linear  averaged  approximation  to  the  dipole  moment 
function,  Eq.  (A3)  in  the  Appendix,  to  facilitate  a  more  parallel  comparison  among  the  various 
theoretical  approaches,  choosing27-29  <mo>— 0.1098  Debye  for  the  averaged  permanent 
dipole  moment  and  <^|>=3.1  Debye/A  for  the  dipole  moment  first  derivative. 
(1  Debye  -  3.336x  10~3"  coulomb-meter.)  The  Newtonian  approach  may  trivially  be  extended 
to  any  dipole  moment  function  and,  at  the  cost  of  some  algebraic  complexity,  the  quantum  and 
7T—0  classical  approaches  may  also  be  extended.27-20-22  Additional  quantum  corrections  to  the 
classical  intensities  may  then  be  required.  Given  a  dipole  moment  function,  any  of  the  theoret¬ 
ical  approaches  can  be  used  to  compute  the  spectral  intensities  for  rotational  and  vibrational- 
rotational  transitions  more  accurately  than  they  are  readily  experimentally  measur¬ 
able.  28.29.33.34 

V.  FROM  GAS  TO  LIQUID 

Having  demonstrated  that  the  quantum-corrected  classical  approach  can  provide  essentially 
correct  spectral  band  contours  for  gas-phase  rotational  and  vibrational-rotational  transitions,  the 
cases  in  which  we  know  accurately  the  potential  energy  and  dipole  moment  functions  and  in 
which  we  can  compare  to  accurate  quantum  calculations,  we  now  turn  to  higher  densities, 
where  our  knowledge  of  the  potential  energy  and  dipole  functions  is  less  certain,  and  our  ability 
to  compute  spectra  by  alternative  means  is  less  developed.  We  will  treat  solutions  of  CO  in 
progressively  higher  densities  of  Ar,  and  finally  in  liquid  Ar. 

Various  approximation  methods,20  classical,  semiclassical  and  quantum,  have  been  applied 
to  computing  vibrational-rotational  band  contours  for  diatomics  in  solution.  Examples  include 
the  stochastic  and  diffusional  approach  of  Bratos  and  co-workers,17-25  the  semiclassical  method 
of  Gordon  and  co-workers,26-28  and  the  classical  impact  study  of  Koszykowski  and  Marcus.20 
Our  choice  to  develop  a  classical  (and  rather  brute-force)  molecular  dynamic  method  which 
includes  all  degrees  of  freedom  is  inspired  by  our  desire  to  extend  the  method  to  more  complex 
solute  and  solvent  molecules  and  to  transient  spectra  in  non-equilibrium  systems. 

Figures  5  and  6  show  experimental  band  spectra  for  room  temperature  CO  in  different 
densities  of  Ar  solvent  after  Coulon  el  a/29  and  in  liquid  Ar  at  a  temperature  of  97  Tv  and  a 
density  of  750  amagats  (or  2.014x  1028  atoms  nT2)  after  Buontempo  et  at.1*  (One  amagat  unit 
of  density  is  the  actual  concentration  of  the  particular  gas  at  OTT  and  1  atm  pressure,  which 
would  be  2.6869xl025  molecules  m-2  for  an  ideal  gas.)40  For  our  theoretical  spectra  we  use  the 
Newtonian  classical  approach  developed  above.  To  reduce  edge  effects  in  our  dynamics  we  use 
minimum  image  periodic  boundary  conditions,  replicating  a  truncated  octahedron  in  a  space¬ 
filling  solid  tessellation,10  smoothly  reducing  our  potentials  to  zero  within  the  radius  of  the 
inscribed  sphere  to  avoid  discontinuities  as  atoms  cross  the  boundaries. 

The  intermolecular  potential  we  use  is  a  simple  pairwise  Lennard-Jones  approximation, 

V(r)  —  4«  {(or/r)12  —  (o-/r)6J,  (37) 

in  which  r  is  the  pairwise  internuclear  distance.  For  Ar  -  Ar  we  use20  a  well-depth  <«124°K  kH 
(equivalent  to  1.712xl0-21  J)  and  a  zero-crossing  radius  <r=0.342  nm.  For  both  Ar  -  C  and  Ar 
-  0  we  use  «— 88 "K  A#  (equivalent  to  I.22x  10~21  J)  and  (r-0.328  nm.  This  Ar  -  CO  potential 
is  a  crude  fit  of  two  identical  two  body  potentials  to  the  potential  surface  of  Parker  and  Pack41 
by  roughly  matching  two  body  potentials  to  the  zero-crossing  and  well  depths  for  the  average  of 
both  three  body  potentials  shown  in  their  Fig.  4.  The  CO  intramolecular  potential  is  as  before. 

In  these  calculations  of  CO  in  Ar,  to  reduce  the  computational  effort  per  specific  spectrum 
and  allow  more  ensemble  averaging,  we  use  in  the  molecular  dynamics  a  step  size  of  10"1'  s, 
and  then  correct  the  spectra  for  the  +16  cm"1  shift  this  produces.  Fifty-five  atoms  of  Ar  are 
used  in  the  unit  cell.  As  can  be  seen  in  Figs.  5  and  6,  the  agreement  of  theory  and  experiment 
all  the  way  from  the  dilute  gas  phase  to  the  liquid  phase  is  quite  satisfactory,  given  the  crude¬ 
ness  of  the  potential  we  have  used,  and  the  experimental  uncertainties. 
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VI.  CONCLUSION 

This  work  shows,  at  least  for  the  simple  system  treated,  that  the  combination  of  classical 
molecular  dynamics,  linear  response  theory,  classical  ensemble  averaging  and  suitable  quantum 
corrections  allows  electric  dipole  rotational  and  vibrational-rotational  band  contours  to  be  com¬ 
puted  which  are  in  close  agreement  with  quantum  calculations  and  with  experimental  results. 
As  we  have  shown,  the  approach  is  equally  applicable  over  the  range  of  densities  from  dilute 
gas  to  liquid  solution  (and  it  should  also  apply  to  solids).  We  have  purposely  avoided,  at  least 
for  now,  the  additional  theoretical  complexity42  inherent  in  treating  the  case  of  strong  absorbers 
in  more  concentrated  solutions,  due  to  the  variation  of  the  index  of  refraction  with  absorption 
frequency. 

The  reason  that  these  classical  calculations  reproduce  so  well  the  observed  spectra  can  be 
understood  as  follows.  For  translational  motion,  the  quantum  numbers  for  the  usual  effective 
volumes  for  atomic  motion  are  large  enough  that  the  correspondence  principle  implies  quantum 
convergence  to  the  classical  limit.  For  rotational  motion,  the  average  quantum  numbers  are 
moderate,  and  if  one  blurs  the  rotational  lines,  for  example  by  collisions,  the  band  contours 
require  at  most  quite  small  quantum  corrections.  For  vibrational  motion,  the  average  quantum 
numbers  are  very  small,  and  the  correspondence  principle  limit  is  far  from  fulfilled.  Thus,  at 
first  glance,  one  is  surprised  that  the  classical  and  quantum  agreement  is  so  good.  The  reason 
for  this  agreement  is  the  unique  equivalence  for  quadratic  potentials  of  many  averaged  measur- 
ables  for  the  classical  and  quantum  cases.42*46  For  example,  the  average  energy  gain  for  the 
forced  harmonic  oscillator  driven  by  a  time  varying  force  Tit)  either  from  radiation  or  from 
collisions  is  the  same  classically  and  quantaily,  and  the  quantum  expectation  values  and  classical 
values  of  position  and  momentum  follow  the  same  classical  equations  of  motion.  Thus  we  find 
exact  agreement  between  our  classical  and  quantum  spectral  calculations  for  harmonic  vibra¬ 
tional  motion  driven  by  the  lime  varying  force  from  the  oscillating  electric  field  of  the  light, 
and  only  the  anharmonic  part  of  the  potential  requires  a  quantum  correction. 

We  are  able  to  ignore  the  effects  of  quantization  of  the  radiation  field,  in  part,  because  the 
lifetimes  for  spontaneous  emission  are  long  with  respect  to  other  changes  in  our  system  once 
we  have  introduced  collisions.  Thus  any  broadening  of  our  spectra  by  radiation  damping  is 
negligible. 

That  classical  mechanics  succeeds  so  well  in  reproducing  the  absorption  spectrum  does  not 
imply  that  the  system  would  behave  classically  if  it  were  examined  by  a  different  set  of  meas¬ 
urements.  Clearly  the  vibrational  properties  are  far  from  classical  if  one  examines,  for  exam¬ 
ple,  the  allowed  energies  of  the  system  instead  of  the  probability  of  light  absorption  versus 
wavelength  or  frequency.  What  is  implied  is  that  we  may  think  about  (or  compute)  the  whole 
system  of  molecules  and  internal  interactions  as  well  as  the  interaction  of  the  system  with  light 
in  a  classical  manner,  applying  our  well-calibrated  and  insightful  classical  intuition,  with  some 
confidence  that  we  can  arrive  at  close  to  the  same  spectral  result  as  if  we  had  properly  handled 
the  real  quanial  nature  of  the  system.  This  opens  the  way  to  understanding  the  vibrational 
spectra  of  much  larger  and  more  complex  systems  under  a  wider  variety  of  conditions  than  is 
feasible  with  a  quantum  approach.  We  thus  believe  that  the  larger  difficulties  in  the  computa¬ 
tion  of  vibrational  spectral  band  contours,  even  for  complex  molecules  and  for  condensed 
phases,  probably  lie  on  the  electronic  side  of  the  Born-Oppenheimer  separation  (potential 
energy,  dipole  moment,  and  polarizability  matrix  as  functions  of  nuclear  position)  rather  than 
on  the  nuclear  motion  side. 

We  plan  in  future  papers  to  show  how  a  parallel  approach  can  be  applied  to  the  electric 
dipole  pure  rotational  and  vibrational-rotational  spectra  of  polyatomic  molecules  as  well  as  to 
Raman  pure  rotational  and  vibrational-rotational  spectra,  both  for  gases  and  condensed  phases. 
In  addition,  we  plan  to  illustrate  the  extension  of  this  molecular  dynamic  and  linear  response 
approach  to  non-equilibrium,  time  dependent  processes,  specifically  to  the  computation  of  tran¬ 
sient  rotational  and  vibrational-rotational  spectra  during  the  course  of  chemical  change.  We 
believe  that  the  computation  and  measurement  of  such  spectra  can  be  a  promising  route  to  dis¬ 
covering  the  molecular  dynamics  of  chemical  reactions  in  solution.2 
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APPENDIX:  QUANTUM  BAND  CONTOURS 

We  will  evaluate  Eq.  (5)  quantum  mechanically  for  a  diatomic  molecule,47  using  rigid- 
rotor  harmonic  oscillator  transition  probabilites,  but  energies  (and  thus  frequencies  and  state 
probabilities)  which  include  anharmonic,  centrifugal  and  vibrational-rotational  coupling  terms. 
The  transition  matrix  elements  are 


</l#a|/>  “  <v'J'M'\p(R)\vJM>.  (Al) 

in  which  v  is  the  vibrational,  7  the  total  rotational  angular  momentum  and  M  the  z  component 
of  angular  momentum  quantum  number  for  a  diatomic  molecule  with  no  spin  or  electronic 
angular  momentum.  Primes  indicate  final  state,  and  we  assume  integration  of  the  dipole 
moment  operator  over  electronic  coordinates  to  give  p.(R),  the  dipole  moment  as  a  function  of 
intemuclear  distance  R.  By  symmetry,  p(R)  must  lie  along  the  internuclear  axis,  so  that 

M(/?) (A2) 

in  which  R  is  a  unit  vector  along  the  intermolecular  axis.  Our  first  approximation  is  to  expand 
/*(/?)  as  a  Taylor’s  series,  keeping  only  the  first  two  terms,,  which  gives  the  standard  linear 
approximation  to  the  dipole  moment  function 


/*(/?)=  no  +  fn\R  (A3) 

in  which  AR  -  (R  -  Re)  is  the  deviation  from  the  position  R,.  at  the  minimum  of  the  poten¬ 
tial  curve  V(R). 

Our  second  approximation  is  to  evaluate  the  matrix  element  in  Eq.  (Al)  using  the  rigid- 
rotor  harmonic  oscillator  approximation47 

< vTA/'l (/*o  +  /*iAR)R|  v7Af>  =  < v'lfyio  +  /*|AR)|v>  <J'M'\R\JM>  (A4) 

-  L>  «v.v  +  ,{(v+l)/2a]4<7A/'|R|7A/>  ,  (A5) 


in  which  8  is  the  Kronecker  delta  and  a  —  m,to,,/ft  (m,  is  the  reduced  mass).  Since  all  M  lev¬ 
els  have  the  same  energy,  we  can  sum  over  M  in  Eq.  (5),  finding48-49 

i  |<7’A/'|R|7Af>|J  -  (7+l)8,_,y  +  J&j+u  .  (A6) 

M—J 


We  can  now  insert  Eqs.  (A5)  and  (A6)  into  Eq.  (5)  and  divide  the  spectrum  into  a  pure  rota¬ 
tional  absorption  spectrum 


/*(<o)  -II  PM)  ^  (7+1)  8 

»-o  J-0 


E(v,J+\)  -  E(v.J) 


ft 


—  u> 


and  a  vibrational-rotational  absorption  spectrum,  in  which  the  R  branch  (7+1—7)  is 


*/,  (co)  -LI  p(v,J)  m,7  (7+1)  I(v+l)/2a)  8 

r»0  J-0 

and  the  P  branch  (7-1—7)  is 


£(v+l,7+l)  -  E{v,J) 


—  cu 


p/i  (io)  -  I  I  p(v,J)  nfj  [(v+I)/2o]  8 

v-0  J- 1 


E(y+1,7-1)  -  E(v.J) 


1i 


—  <o 


(A7) 


(A8) 


(A9) 


where  in  Eqs.  (A7-A9) 


p(v,J) 


_ exp(-/3  E(v,J)\ _ 

f  £,(2y+Dexp(-0£(v,y)) 

v-0  j- o 


(A10) 


The  final  approximation  is  to  use  energy  levels  evaluated  to  second  order  in  perturbation 
theory,  involving  the  third  and  fourth  derivatives  of  the  potential,  giving47 

e(v,j)  -  /t[^(v+i)  +  Beju+ 1)  -  -  a>+i)y(y+n  -  Z>,y2(y+u2J  (aid 


in  which 
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hu} 

4# 


(A12) 

(A13) 

(A14) 

(A15) 

(A16) 

(A17) 


In  the  above,  f'"(7?<,),  )''"'(/?<.)  and  V""(Re)  are  the  second,  third  and  fourth  derivatives  of  the 
internudear  potential  at  the  potential  minimum. 

All  three  of  the  above  approximations  (linearity  of  the  dipole  moment  function,  rigid 
rotor  -  harmonic  oscillator  evaluation  of  the  transition  matrix-elements,  and  energy  level 
evaluation  by  second  order  perturbation  theory  through  third  and  fourth  derivatives  of  the 
potential  function)  may  be  extended  to  higher  terms  to  give  higher  accuracy,33- 50  54  but  in  the 
example  treated  here  it  is  unwarranted. 

Using  these  formulas,  we  can  evaluate  the  quantum  spectrum,  and  then  obtain  the  band 
contours  by  broadening  the  8  functions,  for  example,  into  Gaussians 


8(«yi-w)  —  (■y/n,)Aexp[y(«/,-tu)2]  (A18) 

and  letting  y  decrease  until  the  individual  rotational  peaks  merge. 

By  evaluating  Eqs.  (A7-A9)  at  successively  lower  values  of  Planck's  constant,  until  the 
resulting  functions  converge  to  a  limit,  we  can  compute  the  correspondence  principle  classical 
/f(«),  and  through  applying  the  same  technique  to  Eq.  (2)  we  compute  the  correspondence 
principle  classical  spectrum,  ac  (w). 
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Figure  1.  Absorption  spectrum  for  gas  phase  CO  at  298  TC,  showing  pure  rotational  band  con¬ 
tour  on  left  and  fundamental  vibrational-rotational  band  contour  on  the  right.  Subsequent 
figures  examine  each  of  these  bands  in  detail. 
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Figure  2.  Quantum  corrections  to  classical  vibrational-rotational  spectra.  In  the  upper  panel, 
for  the  real  anharmonic  potential  for  CO,  the  solid  line  (under  the  open  circles)  shows  the  actu¬ 
al  quantum  band  contour  for  gas  phase  CO.  The  dotted  line  shows  the  correspondence  limit 
classical  band.  The  dashed  line  shows  the  effect  of  applying  the  vibrational  anharmonic  ofTset 
quantum  correction  alone  to  the  classical  band  The  dashed  and  dotted  line  shows  the  effect  of 
applying  the  rotational  detailed  balance  quantum  correction  alone.  The  open  circles  show  the 
effect  of  applying  both  corrections  together.  In  the  lower  panel,  similar  contours  are  shown  for 
CO  as  a  harmonic  oscillator,  keeping  the  same  second  derivative  at  the  equilibrium  intcrnuclcar 
distance  as  for  the  real  CO  potential.  There  is  now  no  anharmonic  offset  quantum  correction. 
Note  the  large  increase  in  the  asymmetry  between  the  maximum  intensities  of  the  n  and  R 
branches  produced  by  the  anharmonicity  in  the  potential. 


Figure  3.  Pure  rotational  spectrum  for  gas  phase  CO  at  129K.  The  upper  panel  shows  as  a 
dashed  line  (he  quantum  rotational  lines,  (partially  broadened  by  Eq.  (AI8)  for  illustration)  and 
as  a  solid  line  the  quantum  rotational  band  contour  from  further  broadening  the  individual  rota¬ 
tional  lines  until  they  fuse  into  a  smooth  contour.  The  middle  panel  shows  again  the  quantum 
band  contour  as  a  solid  line  and  the  classical  7)— 0  limit  correspondence  principle  contour  as 
open  circles  which  fall  very  close  to  the  quantum  contour.  The  lower  panel  shows  the  quantum 
contour  as  a  solid  line  and  the  Newtonian  classical  (classical  molecular  dynamics,  classical  linear 
response,  classical  statistical  mechanical  ensemble  average)  band  as  triangles.  The  Newtonian 
spectra  are  an  average  over  an  ensemble  of  10,000  single-molecule  time  histories  of 
24.6xlO~IJs  each  with  I0-,#*s  integration  steps,  run  in  groups  of  300  non-interacting  molecules. 
The  dotted  line  is  the  experimental  gas  phase  spectrum  of  Buontempo  ct  at,  with  Ar  added  to 
partially  broaden  the  rotational  lines.  We  have  multiplied  their  experimental  cross  sections  by  a 
scaling  factor  of  1.15,  which  may  well  be  within  their  experimental  error,  to  better  match  both 
our  calculations  and  the  sum  rule  of  Eq.  (36)  which  they  cite,  both  of  which  arc  based  on  the 
accepted  values  of  the  average  permanent  dipole  moment  of  CO. 
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Figure  4.  Vibrational-rotational  fundamental  band  spectrum  for  gas  phase  CO  at  298  TC.  The 
upper  panel  shows  as  a  dashed  line  the  quantum  vibrational-rotational  lines  (partially  broadened 
by  Eq.  (A18)  for  illustration),  and  as  a  solid  line  the  quantum  rotational  band  contour  from 
broadening  the  individual  vibrational-rotational  lines  further  until  they  fuse  into  a  smooth  con¬ 
tour.  The  middle  panel  shows  again  the  quantum  band  contour  as  a  solid  line  and  the  quantum 
corrected  classical  ft— 0  limit  correspondence  principle  band  as  open  circles  which  match  very 
closely  the  quantum  contour.  The  lower  panel  shows  the  quantum  contour  as  a  solid  line  and 
the  quantum-corrected  Newtonian  classical  (classical  molecular  dynamics,  classical  linear 
response,  classical  statistical  mechanical  ensembte  average)  band  as  triangles.  The  Newtonian 
spectra  are  averaged  over  an  ensemble  of  10,000  single-molecule  time  histories  of  24.6xl0_l-s 
each  with  a  lO'^s  integration  step,  the  dynamics  run  in  groups  of  300  non-interacting 
molecules.  The  dotted  line  is  the  experimental  gas  phase  spectrum  of  Armstrong  and  Welsh 
with  He  added  to  broaden  the  individual  lines,  which  was  given  with  no  absolute  intensity  un¬ 
its,  and  is  thus  scaled  to  match  the  theoretical  curves. 
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Figure  5.  Vibrational-rotational  infrared  fundamental  band  for  CO  in  increasing  densities  of  Ar. 
The  triangles  show  the  theoretical  spectra  for  solutions  of  CO  in  various  densities  of  gas  phase 
Ar  as  computed  from  "Newtonian"  theory  (molecular  dynamics,  linear  response,  statistical 
mechanics  all  carried  out  classically  and  then  quantum  corrected).  Each  spectrum  is  averaged 
over  an  ensemble  of  500  runs  of  24.6x10" l2s  in  length,  with  a  10"l5s  step  size  using  periodic 
boundary  conditions.  The  calculation  of  the  dynamics  for  each  run  took  12  minutes  of  real 
time.  55  Ar  atoms  and  1  CO  are  used,  except  for  the  14  amagat  computation  which  is  1400 
runs  with  17  Ar  atoms.  The  dotted  line  shows  the  corresponding  experimental  spectra  from 
Coulon  c/  al  for  gas-phase  solutions  of  CO  in  Ar  at  room  temperature.  All  the  experimental 
cross-sections  are  scaled  by  a  factor  of  0.82,  which  is  perhaps  wiihin  their  experimental  uncer¬ 
tainty,  and  is  certainly  within  the  range  of  other  historical  experimental  intensity  measurements 
for  the  CO  fundamental  intensity,  in  order  to  more  closely  match  the  theoretical  calculations, 
which  are  based  on  more  recent  dipole  derivative  measurements.  In  all  cases,  the  index  of  re¬ 
fraction,  n,  in  the  linear  response  equations  is  set  to  1.0. 


CO 


Figure  6,  CO  in  liquid  Ar  plotted  to  the  same  scale  as  Fig.  5.  Open  triangles  are  the  “Newtoni¬ 
an*  theoretical  calculations,  for  55  Ar  atoms  and  one  CO  using  periodic  boundary  conditions. 
The  spectra  are  averaged  over  500  runs  of  24.6xl0~12s  each  with  a  10~l:,s  step  size  at  97  TC. 
The  dotted  line  shows  the  liquid  phase  experimental  measurement  of  Buontempo  ei  a I  which  is 
given  with  no  absolute  units,  and  is  thus  scaled  to  match  the  theoretical  calculations.  It  is  ap¬ 
proximately  corrected  for  the  indicated  instrumental  resolution.  The  anhurmonic  shift  comput¬ 
ed  from  Eq.  (35)  is  -25  cm'1  at  97  °K,  which  is  not  quite  sufficient  to  match  the  experimental 
liquid  phase  data,  so  we  have  added  an  additional  -5  cm-:1  solvent  shift  (approximately  a  fre¬ 
quency  scale  factor  of  0.998)  to  our  computed  data  to  match  the  experimental  measurement. 
The  index  of  refraction,  it,  in  the  linear  response  equations  is  set  to  1.0. 
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